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ABSTRACT 

By examining the locations of central black holes in two elliptical galaxies, M 32 and 
M 87, we derive constraints on the violation of the strong equivalence principle for 
purely gravitational objects, i.e. black holes, of less than about two-thirds, rjN < 0.68 
from the gravitational interaction of M 87 with its neighbours in the Virgo cluster. 
Although M 32 appears to be a good candidate for this technique, the high concentra¬ 
tion of stars near its centre substantially weakens the constraints. On the other hand, 
if a central black hole is found in NGC 205 or one of the other satellite ellipticals of 
M31, substantially better constraints could be obtained. In all cases the constraints 
could improve dramatically with better astrometry. 


1 INTRODUCTION 


The strong equivalence principle (SEP) states that any ob¬ 
ject regardless of its composition will travel through a grav¬ 
itational field in the same way. This includes objects with 
varying contributions of gravitational energy. In particular a 
black hole whose mass is entirely gravitational should travel 
in the same manner through the gravitational field as a star. 
[Nordtve'dt] ( 1968a|b I argued that metric theories of gravity 
other than general relativity may exhibit violations of the 
SEP — in particular, objects which have a large contribu¬ 
tion of gravitational energy to their makeup may travel dif¬ 
ferently through a gravitational field than other objects. 

Typically the strong equivalence principle is probed by 
looking for the Nordvedt effect in the Earth-Moon system or 
binary stars consisting of a neutron star and a white dwarf 
( Stairs et al.||20d5 K The Nordtvedt (1968cI effect results in 
the polarization of an orbit in the direction of the gravi¬ 
tational acceleration of a large body: the Sun in the case 


of the Earth-Moon system (Williams et al. 1976 Shapiro 


|et al.|[l976p , and the Galaxy in the case binary pulsars. 
The recently discovered millisecond pulsar in a triple sys¬ 
tem, PSR J0337-I-1715, will provide further interesting con¬ 
straints ( [Ransom et al.|2014[ ). 

Here we perform the test proposed by |Hui fc Nicolisj 


(20121. In particular we will examine the polarization of the 


orbits of supermassive black holes through the central re¬ 
gions of elliptical galaxies. We will focus on elliptical galax¬ 
ies because the presence of central black holes is common in 
sufficiently large elliptical galaxies and the location of the 
bottom of the gravitational potential is straightforward to 
constrain by determining the centroid of the distribution of 
stellar light. Interactions between galaxies are ubiquitous, so 
example systems are straightforward to find. In particular 
we will examine the small elliptical galaxies that orbit the 
largest neighbour of the Milky Way galaxy, the Andromeda 
galaxy or M31, and the massive elliptical galaxy, M87, in 
the Virgo cluster. We will derive an upper limit on SEP vi¬ 


olation - how this limit should be interpreted in the context 
of general gravitational theories will be discussed in © 


2 CALIBRATION 


Let us suppose that we are studying a small stellar system in 
the gravitational field of a large one. The stars in the small 
system experience an acceleration toward the centre of the 
larger one 


GM 


a* 




( 1 ) 


where d is the distance between the small system and the 
larger one and M is the mass of the larger system. These ac¬ 
celerations are typically ten to one hundred times larger than 
those exerted by large-scale structure a ~ 600 km s~^Ho ~ 
10"^°cm s“^ 

Furthermore, let us assume that a massive black hole 
within the smaller system experiences a different accelera¬ 
tion 


a. = (1 — A) a* 


( 2 ) 


where A quantifies the extent of violation of the SEP. 


Nordtvedt (19821 presents an equivalent definition where the 


ratio of the inertia masses and passive gravitational mass of 
an object could depend on its constitution and in particular 
the contribution of gravitational energy. 


m-Gp Ua 

- = 1 - riN -j. 

mi mc^ 


(3) 


In the context of the parametrized post-Newtonian treat¬ 
ment of deviations from general relativity ( Will|2014 1, 

10 2 2 1 

VN = 4:l3-^-3-—£,-ai +-a2--Cl--(2- (4) 


For the case of a black hole, Uq = me? so A = rjN- 

The massive black hole also experiences an acceleration 
due to the stellar system in which it resides. Furthermore, 


© 2015 RAS 





























2 Asvathaman, Hey I & Hui 


dynamical friction with the less massive stars nearby damps 
any motion of the black hole relative to its equilibrium po¬ 
sition on a short time scale, so absent SEP violation or rare 
astrophysical effects such as a recent merger or asymmetric 
jets the black hole should lie at the bottom of the gravi¬ 
tational potential well. The effects of individual stellar en¬ 
counters, Brownian motion, are much smaller than what we 
discuss here ( Broderick et al.|2011 l. 

Near the minimum of the stellar potential we can ap¬ 
proximate it as a harmonic oscillator so 


NGC 147 ■ 


(1 — A) oj'^x = a, — a-i, = —A 


GM 


( 5 ) 


and X is the displacement of the black hole from the centre of 
the small stellar system in the direction of the larger system. 
If A > 1 the black hole would not lie near the centre of the 
galaxy, so this region is excluded. 

The value of ui depends on the structure of the small 
stellar system and depends on the mean density within the 


region surrounding the black hole. For a Hernquist (19901 


model for the density distribution that well describes the 
stellar distribution of elliptical galaxies 


P = 


1 


27r r (r -I- a)^ 


( 6 ) 


where RbS, the half-light radius in projection is approxi¬ 
mately 1.8153a. We take the black hole mass to be m, to be 
much less than that of the small system. However, within 
a small region near the black hole, the mass of the black 
hole dominates that of the stars. We will take the mass of 
the stars within this sphere of influence to be am, where 
a « 1. Now if the black hole is displaced from the centre 
by violation of the equivalence principle, Eq. then we as¬ 
sume that these stars will be displaced with the black hole 
and remained centred on the black hole with a Hernquist 
distribution of density, Eq. On the other hand, we will as¬ 
sume that the rest of the galaxy will remained centred on the 
undisplaced position of the black hole. Of course, there will 
be a smooth transition between these two regimes, but this 
approximation contains the essential point that the cusp of 
stars at the centre of the galaxy will remain centred on the 
black hole. By numerically integrating over the gravitational 
potential energy of this configuration we find 


; 0.30 

V am, J 


1/2 Gm 


( 7 ) 


where the coefficient comes from the numerical integration 
and the scaling comes from considering the structure of the 
Hernquist model at small radii. 

If one directly used Eq. © to calculate the restoring 
force toward the centre, one would find that it is constant 
with displacement from the centre. This would assume that 
the density is actually singular at the centre but also would 
neglect the effect of the black hole on the neighbouring stars 
which would be dragged with it. We will compare the esti¬ 
mate including the effect of neighbouring stars (Eq. with 
the detailed measurements of the potential near the core of 
the elliptical galaxy, M 87, as inferred from stellar velocities 
and find good agreement. 

We will take a = 1 and m, = 0.005m where m is the 
mass of a bulge and halo approximately modelled by a Hern¬ 
quist model (as in|van der Marel et al.]1997b|, we find that 


NGC 185' 


M32/ /M31 


■NGC ^05 


Figure 1. The geometry of the largest elliptical satellite galaxies 
of the Andromeda Galaxy. The Sun is to 773 kpc to the left of 
M32, and the radius of the dashed circle is 100 kpc. 


^=0.24-A^^-, (8) 

d 1 — A m 

so the relative displacement only depends on the relative 
masses and sizes of the various bodies and the angular dis¬ 
placement does not depend on the distance to the systems 
from Earth. We can make this explicit by calculating the 
apparent displacement across the sky as 


Ae = 0.24 


A M a^ 

1 — A 52 ui 


3 . 
cos ^ 


( 9 ) 


where i is the inclination of the line connecting the two bod¬ 
ies with respect to the plane of the sky, 5 is the angular 
distance between the bodies and a is the angular size of the 
Hernquist radius a of the host galaxy. If smaller galaxies fol¬ 
low an isothermal distribution about the perturbing galaxy 
as Koch & Grebel ( 2006| found for the satellites of M 31, the 
mean value of the geometric term | cos^ i\ is 4/(37r) ~ 0.42 
and the median value is \/2/4 « 0.35. Ninety percent of the 
time the term lies between 5 x 10“^ and 0.99, so geometric 
information is especially useful to give firm constraints. Of 
course the halo of the galaxy group must end somewhere so 
the lower limit on cos® i found here is somewhat unrealis¬ 
tic. On the other hand if the density of satellites is propor¬ 
tional to r“® like the outer regions of the Navarro-Frenk- 


White proHle (19961, the mean value of the geometric term 
is 37r/16 « 0.59 and the median value is 3-\/3/8 « 0.65. 
Ninety percent of the time the term lies between 0.03 and 
0.996. 


3 THE ANDROMEDA GALAXY SYSTEM 

The nearest neighbour of the Milky Way has four small el¬ 
liptical galaxies orbiting it: M32, NGC 205, NGC 147 and 
NGC 185. Fig. [T] depicts the location of the various galaxies 
relative to each other from the distances compiled by |Koch] 
& Grebel ( 2006[ ). Of course, the distances to the various 
galaxies are more uncertain than their positions on the sky. 


A supermassive black hole (SMBH) has already been 
discovered in the core of the compact elliptical galaxy M 32 
(van der Marel et al. 1997aI and found to be within 0.5 
arcseconds of the nucleus of the galaxy (Yang et al.||2015|. 
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M32 is located about 0.4° from the centre of M31 in the 
plane of the sky. We can estimate the displacement of the 
black hole from the centre of M 32 as a function of A in the 
best-case scenario of the minimal separation of M 32 from 
M 31 of about 5 kpc. In this case M 32 would be well within 
the halo of M31, and we could use the circular velocity of 
M31 of 200 km/s at this distance to estimate the acceleration 
of M32 toward the centre of M31, a* in Eq. Q, 

2 

a* = ^ = 3.6 X 10~®cm s~^. (10) 

For the galaxy M 32 as a first approximation we will use 
a Hernquist model that fits to the global structure of the 
galaxy. This model fits the observed half-light radius of the 
galaxy (we will revisit this assumption later in light of data 
that probe the centralmost regions, see Fig. [^, so we take 
m = 8 X 1O®M0 and Retf — 100 pc so a = 55 pc and 


X = 0.091—^yrpc 

1 - A 

(11) 

or 


A9 = -— - = 0.0288 ^ arcseconds, 

aM32-Sun 1 — A 

(12) 


so if M 32 lies indeed at the close to minimum distance to 
M31, the existing measurement of the position of the black 
hole yields a constraint of A < 0.95. 

We can improve upon the astrometric constraint of 
[Yang et al.| ( |2015| ) if we assume that the position of the black 
hole is coincident with the cusp in the surface brightness 
of the galaxy. This is a reasonable assumption because the 
black hole would entrain the neighbouring stars, so the stel¬ 
lar nucleus should be centred on the black hole. To measure 
the possible deviation of the position of the black hole on 
the sky and the centre of the potential well of the galaxy we 
use images from the Wide-Field Camera 3 (WFC3) on the 
Hubble Space Telescope (HST) in the filter, F555W. This 
instrument features 0.04 arcsecond pixels and a 162 x 162- 
arcsecond field of view. Both the Space Telescope Imaging 
Spectrograph (STIS) and the planetary camera on Wide- 
Field Planetary Camera 2 (WF/PC2) feature a finer pixel 
scale but at the expense of field of view. We used images from 
GO-11714 (PI: Bond) with total exposure time of 120 sec¬ 
onds. To localise the stellar core which lies within the black 


hole sphere of influence of about 0.3 arcseconds {Joseph et al. 


2001), we used circular and elliptical apertures of 0.16 arc- 


seconds and 0.2 arcseconds. To localise the centroid of the 
outer regions we used two elliptical annuli: one of inner axes 
2.8 and 4.3 arcseconds and outer axes of 5.7 and 8.5 arc- 
seconds and one half that size, well outside the black hole’s 
sphere of influence. The centres of the annuli are initially 
roughly placed within about an arcsecond of the core and 
updated to lie on the calculated centroid. This is iterated 
thirty times. The centroids of the various annuli to probe the 
galactic potential coincide to within 0.01 arcseconds as do 
the ellipses and circles centred on the stellar cusp, yielding 
an estimate of the precision of 0.01 arcseconds. Furthermore, 
the centroid of the cusp and the outer regions also agree to 
within 0.01 arcseconds. If we combine this constraint with 
the distance estimate of Koch & Grebel (20061, this yields a 
constraint —0.53 < A < 0.26. In fact the maximal distance 
between M 31 and M 32 consistent with the uncertainties 


given by Koch & Grebel (20061 is 45 kpc (see also Dierickx 



Semimajor Axis [arcseconds] 


Figure 2. The surface brightness in the core of M32 from |Lauer| 
|et al.| | |1998[ ) , [Michard &: Niet^ ( |1991^ and the photometry and 
models used in this paper. We assumed B — F555W = 1 to place 
the|Michard &: Nieto]|[l991[) data on this plot. 


et al.|20l4 for an alternative interpretation). In the case the 
constraints on A would be much weaker with A < 0.99. 

Before we focus on the smaller elliptical galaxies of the 
M31 system, we will examine the central region of M32 in 
further detail. |Lauer et al.| (|1998|) deconvolved images with 
HST and WFPG-1 and WFPG-2 and found a very high cen¬ 
tral concentration of stars. Fig. [^depicts the surface bright¬ 
ness in the core of M 32 from ground-b ased data ([Michard &:| 
Nieto|19^ I and HST WFPG2 images ( [Lauer et al. 1998) as 
well as the WFG3 images and the models that we are using in 
this paper. We can see that the global Hernquist model that 
reflects the half-light radius of the galaxy of about 27 arc- 
seconds dramatically underpredicts the light in the core of 
the galaxy. On the other if we focus on the inner two arc- 
seconds of the galaxy we can characterise the distribution of 
light with a much more concentrated Hernquist model with 
a = 1.3 arcseconds or 5.5 pc. The total mass in this cen¬ 
tral Hernquist model is about 0.17 of the entire galaxy or 
1.4 X 1O®M0. This more concentrated model underpredicts 
the light outside of two arcseconds, but it does well in the 
region of the central black hole whose mass we have assumed 
to be about 4 x lO^M©, still significantly smaller than the 
mass of the central Hernquist model, so we can apply the 
preceding analysis here but with the new parameters for the 
Hernquist model. 

If we examine Eq. we see that the spring constant 
of the Hernquist model that best reflects the central region 
of M 32 is about 70 times larger than that of the model of 
the entire galaxy and the expected displacement is 70 times 
smaller, 


A6» = 0.4 


1 - A 


milliarcseconds. 


(13) 
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The current astrometry yields A6 < 10 milliarcseconds, 
yielding a weaker constraint of A < 0.96. 

The three smaller elliptical galaxies, NGC 205, NGC 147 
and NGC 185, have smaller mass densities and may have 
more favourable geometry. However, SMBHs have not yet 
been identihed in these galaxies. Each of these galaxies has 
a mass of about 2 x IO^Mq and a i?eft ~ 300 pc, so a ~ 
165 pc. The two galaxies, NGC 147 and NGC 185 are each 
about 100 kpc from M31, so the enclosed halo mass is about 
IO^^Mq. This yields a typical value of the displacement of 

X = 0.55—^-^pc and A9 = 0.18—arcseconds.. (14) 

Although NGC 205 may be somewhat closer to M 31 at 
50 kpc than the other galaxies, the geometry is somewhat 
poorer according to the distance given by |Koch fc Grebel| 
120061 perhaps reducing the angular displacement by a fac¬ 


tor of two counteracting the increased acceleration. This 
yields a similar observed displacement. If a SMBH is dis¬ 
covered near the centre of these galaxies and can be lo¬ 
calised within 0.5 arcseconds, this would yield a constraint 
of A < 0.74. 

On the other hand, NGC 205 may lie at the same dis¬ 
tance from us as M31. This is consistent with the errorbars 
given by [Koch fc Grebel| ( |2006| ). Furthermore, [Geha et ah] 
(20061 argue that NGC 205 is interacting tidally with M31, 


favouring this interpretation. At this minimum possible dis¬ 
tance to M 31 of about 7 kpc, and we obtain a much larger 
linear and angular displacement 


A ^ A 

X = 10.:^-:-pc and A0 = 3.1--arcseconds. 


'1 - A' 


A - A 


(15) 


In this most favourable case the constraint would be —0.19 < 
A < 0.14 with 0.5-arcsecond astromety and |A| < 3 x 10“^ 
with 0.01-arcsecond astrometry. 


4 VIRGO ELLIPTICAL GALAXIES 

One of the largest elliptical galaxies in the Virgo Cluster, 
M 87, also contains a supermassive black hole and lies at 
the centre of the Virgo A subcluster which contains the bulk 
of the mass of the Virgo Cluster, about IO^^Mq. The Virgo 
Cluster is still in the process of forming and contains at least 
two additional subclumps of masses of about lO^^M©. The 
subclump closest to M 87 is centred on M 84 and M 86 about 
1.3° away. M89 is also a similar angular distance away, so 
this analysis applies to it as well, but it is likely to be less 
massive. 


For the galaxy M87 we take R^s = 164 (Ferrarese 
et al.|2006 1 so a = 90”. We will take the mass of M 87 to be 
about 4 X IO^^Mq | Wu fc Tremaine|2006 1 and take the mass 
of the subclump containing M 84 and M 86 to be lO^^M© at 
a distance of 1.3°. This yields a displacement of 


A0 = 0.02 cos° i 


1 - A 


arcseconds. 


(16) 


We can examine the stellar distribution of M87 in fur¬ 
ther detail. We use images from the Wide-Field Camera 3 
(WFC3) on the Hubble Space Telescope (HST) in the ultra¬ 
violet filter, F275W. We used images from GO-12989 (PI: 
Renzini) with total exposure time of 5599 seconds as de¬ 
picted in Fig. 1^ M 87 is an EO galaxy so its isophotes are 



Figure 3. Central region of M 87 focusing on the AGN. Several 
apertures to determine the centroid of the light are depicted. The 
two middle annuli lie from 9 to 15 arcseconds and 12 to 24 arcsec¬ 
onds from the centre. The inner apertures are circular and 0.15 
and 0.2 arcseconds in radius. The outermost aperture between 24 
and 36 arcseconds is not depicted. The galaxy M 89 is east (left) 
of M 87 and M 84 is to the West (right). 


nearly circular ( Ferrarese et al.|2006 1, so to measure the sur¬ 
face brightness we average over circular annuli of width of 
one arcsecond that exclude the jet and the direction oppo¬ 
site the jet as depicted by the regions in Fig. (but with 
larger widths). Fig. [^depicts the resulting surface bright¬ 
ness profile that extends to 70 arcseconds, the edge of the 
image. This is well within the half-light radius of the galaxy, 
so we cannot probe the distant light well, and the surface 
brightness measured at large radii depends on the assumed 
sky brightness which we also ht in parallel with the Hern- 
quist model. We again see that the global Hernquist model 
based on the half-light radius does not fit the central region 


well. Ferrarese et al. (20061 find that the central region of 
M87 is better fit by a “core-Sersic” model (Sersic|19M 


Tru- 


jillo et ah||2004l that has a much less cuspy central surface 


brightness distribution than a Hernquist model. The break 
radius between the core and the outer Sersic model is about 
7 arcseconds. Furthermore, in the core of M87 we have de¬ 
tailed kinematic data, so we can estimate the properties of 
the gravitational potential induced by the stars. 


Walsh et al. (20131 determine the mass of the central 
black hole of M87 to be about 3.5 x 1O®M0 (approximately 
10“® of the mass of the galaxy) and to dominate the mass 
within about 5.61 arcseconds of the centre. They also model 
the contribution of the stellar mass to the circular velocity 
in this inner region, yielding a value of about 160 km/s at 4" 
from the black hole consistent with the photometry of |Fer-| 
rarese et al. (20061 and a mass-to-light ratio of 4. Assuming 


a smaller mass-to-light ratio would yield smaller constraints 
on A, roughly the expected displacement for a given value of 
A is inversely proportional to the square of the mass-to-light 
ratio. As the stellar mass-to-light ratio decreases, the effect 
of the stellar potential diminishes, increasing the displace¬ 
ment, and the sphere of influence of the black hole increases, 
so we can probe the effect at larger radii from the black hole, 
further increasing the expected displacement. 

Within 10 arcseconds the black hole contributes about 
thirty percent of the mass, so we will use the circular ve¬ 
locity due to the stellar contribution at this radius of about 
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Radius [Arcseconds] 


Figure 4. The surface brightness in the core of M 87 as measured 
from the image depicted in Fig.[^along with the Hernquist model 
based on the measured of 164” and one based on the inner 
F275W photometry with R^g of 74". 


210 km/s to estimate the displacement of the black hole due 
to SEP violation. Using a distance to M 87 of 16.4 Mpc (Bird 
|et al.|2010p and the same mass and angular distance to the 
subclump centred on M 84 as earlier yields 

A 


X = 1.1 


and 


1 - A 


cos I pc 


(17) 


(18) 


A 3 

A0 = 0.014-- cos i arcseconds 

1 - A 

where we have extrapolated the theoretical stellar mass pro¬ 


file of Walsh et al. (20131 out to 10 arcseconds, slightly 


beyond the break radius of 7.15 arcseconds determined by 


Ferrarese et al. (20061 where the surface brightness profile 
steepens. 

The position of the black hole can in principle be de¬ 
termined to microarcsecond precision with microwave inter¬ 


ferometry (Broderick et al. 20111. On the other hand, the 


centre of the potential well is best defined in the optical 
using the isophotes of the galaxy. Fortunately, the super- 
massive black hole is also apparent in the optical. The in¬ 
nermost isophotes will be centred on the black hole because 
it will dominate the potential well, so the centres of more 
distant isophotes provide an estimate of the potential well 
that constrains the black hole. 

Batcheldor et al. (2010| argued from isophotal analysis 


of observations with the Advanced Camera for Survey on 
HST that the AGN may be displaced from the centre of the 
galaxy by about 7 pc in a direction opposite to the observed 
jet (about 0.1 arcseconds). However, subsequently neither 


Gebhardt et al. (20111 nor Walsh et al. (20131 found evl 


dence for displacement greater than about one parsec; the 
position determined by Gebhardt et al. (20111 was consis¬ 


tent with that of Batcheldor et al. (20101 within the larger 


errorbars of Gebhardt et al. (2011 Batcheldor et al. (20101 
outlined several astrophysical explanations for a potential 
displacement such as a SMBH binary, a recent merging of 
black holes or a one-sided jet. The interaction with the 
neighbouring stars or even clusters of stars is too weak to 
explain such a displacement. 

To measure the possible deviation of the position of the 
black hole on the sky and the centre of the potential well of 
the galaxy we calculated the light centroid within three cir¬ 
cular annuli of inner and outer radii 24 and 36 arcseconds, 12 
and 24 arcseconds, and 9 and 15 arcseconds. If we take the 
mass of the black hole to be about 4 x 1O®M0, it will domi¬ 
nate the mass within about 6.6 arcseconds from the centre. 
Within 12 arcseconds it contributes less than 25% of the 
mass, and within 24 arcseconds it contributes about 8% of 
the mass; therefore, these annuli lie outside the black hole’s 
sphere of gravitational influence. The slices of the annulus 
along the direction of the observed jet and potential coun¬ 
terjet are omitted from the calculation and so exclude light 
from the AGN itself. This technique contrasts with that of 


Batcheldor et al. (20101. They explicitly masked the notice¬ 


able jet emission as well as globular clusters. Here we remove 
a much larger region from the analysis both in the direction 
of the jet and the opposite direction. This excludes both 
the observed jet and a wide regions around the jets. The 
excluded region is symmetric, so we minimize any potential 
bias along the jet axis. Furthermore, the measurements that 
we use are in F275W with a higher angular resolution cam¬ 
era (WFC3 vs. ACS). In this ultraviolet filter, the emission 
from globular clusters is negligible. The centre of the annuli 
are initially roughly placed within about an arcsecond of the 
AGN and updated to lie on the calculated centroid. This is 
iterated thirty times. 

To determine the location of the black hole we per¬ 
formed a similar centroiding on circular apertures of 0.15 
and 0.2 arcseconds, also initially centred on the AGN within 
0.2 arcseconds. We have repeated this process several times 
using different starting positions, resulting in slightly dif¬ 
ferent centroids all consistent within their mutual standard 
deviation of 0.03 arcseconds, slightly less than a pixel. This 
is consistent with the precision found by Walsh et al. (20131 
and poorer than we found with M 32 because the surface 
brightness profile of M 87 is much shallower. The isophotes 
constrain the deviation of the black hole from the centre of 
the potential within 0.03 arcseconds, yielding a constraint on 
the violation of the SEP for black holes of A < 0.68, assum¬ 
ing the more conservative displacement estimate Eq. (181 
and a favourable geometry of i « 0°. Using additional im¬ 
ages, more sophisticated isophote fitting procedure and more 
detailed modelling of the mass distribution of the Virgo clus¬ 
ter could possibly yield smaller constraints. 


5 CONCLUSIONS 

How does our limit on A fit within the context of existing 
scalar-tensor theories? SEP violation occurs in these theories 
because compact objects have a suppressed scalar charge, 
leading to a weaker scalar coupling to the environment. For 
the classic Brans-Dicke theory, solar system tests constrain 
the scalar-matter coupling to be fairly weak to begin with 
(much weaker than tensor-matter), thus predicting SEP vi- 
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olation that is much below our limit. More recent scalar- 
tensor theories open up the interesting possibility of predict¬ 
ing observable SEP violation on large scales while respecting 
the solar system constraints - the latter is achieved by what 
is known as screening mechanism. Screening introduces a 
subtlety in the interpretation of A, however - it becomes 
scale/environment dependent. In the case of the M 32-M31 
system, we estimate A ~ 10~^ in galileon theories (Nicolis 


et al. poo^ . For snch theories, A ~ 1 can be found in situa¬ 


tions where the acceleration is due to large scale structure. 


quite a bit smaller than the acceleration of our setup (Hui 
fc Nicolis||2012 1. Despite this, deriving a phenomenological 


limit on A is useful, for several reasons]^ First, new theories 
- massive gravity being an example - are still being devel¬ 
oped. Second, equivalence principle violation remains to be 
worked out even for some existing theories ( Hui et al.|2009 1. 
Third, this type of galaxy-satellite system constraint on A 
can be improved. 

Both the calibration and observational analysis pre¬ 
sented here can be improved upon in several directions. From 
the calibration standpoint, we have considered just the force 
on the black hole and the coterie of stars in its immediate 
gravitational influence. Moving the black hole from the cen¬ 
tre of the stellar system will also move these stars, and we 
have used a simple model for their coupled motion. Addi¬ 
tional study is warranted where one would start the black 
hole in the centre of the galaxy increase the differential ac¬ 
celeration gradually, say due to the effect of the small galaxy 
approaching the larger one, and determine where the black 
hole ends up. 

As we have mentioned earlier, there are more sophis¬ 
ticated ways to determine the centre of the galaxy with 
isophotal htting and by considering measurements from sev¬ 
eral observations at several wavelengths. In principle one 
could find tighter constraints on the deviation of the black 
hole from the centre of the potential. In the case of M 32 the 
existing constraints are just 0.5 arcseconds because one has 
to link the radio coordinate system to the visual one. These 
could possibly be improved dramatically with further obser¬ 
vations. The visual core of M 32 lies within 0.01 arcseconds of 
the centroid of the outer isophotes of the galaxy, indicating 
that if the black hole is coincident with the core, rjN < 0.96, 
using a density model adapted to ht the central regions of 
M32. The discovery of a supermassive black hole in one of 
the other elliptical satellites of M31 could provide stronger 
constraints from the M 31 system. However, the main uncer¬ 
tainty in the constraints comes from the unknown geometry 
of the system and the central potential as we demonstrate 
for M 32 in particular so further information on the relative 
distances and the central kinematics of the various galaxies 
would yield better constraints on the results. Observations 
of the black hole in M87, the stellar kinematics near the 
black hole and an assumption of a favourable geometry for 
the neighbouring galaxies within the Virgo cluster yields a 
constraint of rjN < 0.68, competitive with the results from 
pulsar timing ( Stairs et al.|2005 l. Perhaps other interacting 
and nearly interacting galaxies in the local Universe could 


^ Note also the environmental dependence of A is often such that 
A on the left hand side of Eq. © is much smaller than A on the 
right; the statements above refer to the latter. 


also provide more information and would statistically limit 
the violation of SEP in the motion of black holes. 

The constraint on SEP from the motion of black holes 
probes the SEP in an essentially different limit from that of 
lunar ranging; however, the technique outlined here could 
become competitive with the lunar ranging results which 


now constrain |7;jv| < 1.3 x 10 ® (Baefiler et al.|[r999 Adel- 


berger 20011. Improved astrometric measurements of the 


centre of the potential of M 87 and the location of the black 
hole could be achieved with existing visual and UV data 
and the position of the black hole could be constrained fur¬ 
ther with radio interferometry provided the optical and radio 
astrometry could be tied together with sufficient precision. 
Relative astrometry with a precision ~ 10“'* arcseconds can 
be achieved with Hubble photometry (e.g. Heyl et al.|2012 1, 
a factor of one hundred better than presented here, so con¬ 
straints on |?7iv| of 10“® could be possible with supermassive 
black holes. If microarcsecond astrometry can be brought to 
bear (e.g. [Broderick et al. 2011[ ) on M87, then constraints 
on r]N on the order of 10~^ might be possible. 
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